The Astrophysical Journal, scheduled to vol. 628, no. 2, 2005 July 20 

Preprint typeset using 1^1^^ style cniulatcapj v. 4/12/04 



in 

o 
o 

(N 



> 

o 
o 
\o 
o 

o 



Oh 
I 

o 



X 



RADIATION FROM CONDENSED SURFACE OF MAGNETIC NEUTRON STARS 
Matthew van Adelsberg\ Dong Lai\ Alexander Y. Potekhin^"*, Phil Arras"* 

Received 2004 May 28; accepted 2005 April 11 

ABSTRACT 

Recent observations show that the thermal X-ray spectra of many isolated neutron stars are fea- 
tureless and in some cases (e.g., RX J1856.5— 3754) well fit by a blackbody. Such a perfect blackbody 
spectrum is puzzling since radiative transport through typical neutron star atmospheres causes no- 
ticeable deviation from blackbody. Previous studies have shown that in a strong magnetic field, the 
outermost layer of the neutron star may be in a condensed solid or liquid form because of the greatly 
enhanced cohesive energy of the condensed matter. The critical temperature of condensation increases 
with the magnetic field strength, and can be as high as 10^ K (for Fe surface at i? ^ 10^"^ G or H 
surface a.t B ^ a few x lO^"* G). Thus the thermal radiation can directly emerge from the degenerate 
metallic condensed surface, without going through a gaseous atmosphere. Here we calculate the emis- 
sion properties (spectrum and polarization) of the condensed Fe and H surfaces of magnetic neutron 
stars in the regimes where such condensation may be possible. For a smooth condensed surface, the 
overall emission is reduced from the blackbody by less than a factor of 2. The spectrum exhibits 
modest deviation from blackbody across a wide energy range, and shows mild absorption features as- 
sociated with the ion cyclotron frequency and the electron plasma frequency in the condensed matter. 
The roughness of the solid condensate (in the Fe case) tends to decrease the reflectivity of the surface, 
and make the emission spectrum even closer to blackbody. We discuss the implications of our results 
for observations of dim, isolated neutron stars and magnetars. 

Subject headings: stars: magnetic fields — radiation mechanisms: thermal — stars: neutron — X-rays: 
stars 



1. INTRODUCTION 

In the last few years, much progress has been made in 
studying surface radiation from isolate d neutron stars 
(NSs) (see, e.g.. iPavlov fc Zavlinll2?M for a review). 
So far about 20 NSs have been detected in thermal 
emission. With the exception of 3-4 sources,^ the 
thermal spectra of most observed isolated NSs are fea- 
tureless and sometimes well fit by a blackbody. For 
example, deep observations with Chandra and XMM- 
Newton show that the soft X-rav (0.15-1 key ) spec- 
trum of RX J1856. 5-3754 l|Walter et aUlTOQ^ can be 
fit wi th an almost perfect blackbody at kT = 64 eV 
(e.g. , IDrake et a,].l Hool IBurwitz et alJ ^QQ^ . The op- 
tical data of RX J1856.5— 3754 is well represented by 
a Rayleigh-Jeans spectrum, but the observed flux is a 
factor of 7 hig her than ex trapolation from the X-ray 
blackbody (see iPons et aLl 120021) . Thus the spectrum 
of RX J1856.5— 3754 is best fit by a two-temperature 
blackbody model. Using this model as well as the ob- 
servational upper limit (1.3% at 2(7) of X-ray pulsa- 
tion (Burwitz et al. 2003), Braie & Romani (2002) ob- 
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tained several constraints on the viewing geometry, mass- 
to-radius ratio, and temperature distribution. Another 
much-studied, dim, isolated NS, RX J0720. 4-3125 also 
shows an X- ray spectrum well fit by a blackbody at 
T - 1 MK l|Paerelset al J 1200 It but see IHaberl et alJ 
l2004al for possible spectral features). 

The featureless, and in some cases "perfect" blackbody 
spectra observed in isolated NSs are puzzling. This is 
because a NS atmosphere, like any stellar atmosphere, is 
not a perfect blackbody emitter due to nongrey opaci- 
ties: On the one hand, a heavy-element (e.g., Fe) atmo- 
sphere would produce many spectral lines in the X-ra; 
band (e.g., Raiagooal & Romani 1996; Pons et al. 200 
on the other hand, a light-element (e.g., H or He) atmo- 
sphere would result in an appreciable hard tail relative 
to the blackbody fe-g.- lShibanov et al]ll992D . 

One physical effect that may help explain the ob- 
servations is vacuum polarization. Recent work has 
shown that for surface magnetic fields B > 10^** G, 
strong-field vacuum polarization can significantly affect 
radiative transfer in NS atmospheres, leading to depres- 
sion of the hard spectral tail and suppression of the 
(cyclotron or atomic) absorption lines ( Lai & Ho 200^ 
'2003a; Ho & Lai 2003, 2004: 'IIolt aL 2003; Llovd 2003). 
Indeed, Ho & Lai (2003 1 suggested that the absence of 
lines in the observed the rmal spectra of s e veral anoma- 
lous X-ray pulsars (e.g . . iJuett et al.|l200a iTiengo et alJ 
2002; Mor net alJl2003t iPatel et alJl2003F can be natu- 
rally explained by the vacuum polarization effect. For the 
dim isolated NSs RX J1 308.6-f2127 (Haberl et alJl2003D 
and RX J1605.3-F3249 ijvan Kerkwiik et al.ll2004|) . the 
observed line features are consistent with surface fields 
< 10^"* G, at which vac uum polarization does not affect 
the emergent spectrum (|Ho fc Laill2004j) . In the case of 
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RX J1856.5— 3754, if t he NS has a magnetar-l ike sur- 
face magnetic field fsee lMori fc R.udermanll2fl()3 ^ it may 
be possible to explain the almost perfect X-ray black- 
body with an atmosphere. However, theoretical mod- 
els for low-temperature (kT ~ 60 eV) magnetar atmo- 
spheres are currently not available because of uncertain- 
ties in treating atomic, molecular opacities a nd dense 
plasma effects in such cool atmospheres (see iHo et all 
|S)03; Potckhi n'&"Chabrieii r2003 . 200|). 

Recently, several groups have suggested that the spec- 
trum of RX J1856.5— 3754 might be explained if the 
NS has a condensed su r face with no atmosphere above 
it llBurwitz et alJ IpM l2003t iMori fc RudermanI l20nl 
iTuroUa et alJl200^ . The notion that an isolated mag- 
netic NS has a condensed surface was first put forward 
in the 1970s fsee Rudcrman 1971: Flowers et al. 1977| 
although these early studies overestimated the cohesive 
energy of Fe solid at B ~ 10^^ G. Revise d cal culations 
yielded a much smaller cohesive energy (|Miiller 1984 
|jlpnes"198(ij iNeuhauser et al .''1987^ making condensation 
unlikely for most observed NSs. Lai & Salpetcr ( 1997) 
studied the phase diagram of the H surface layer of a NS 
and showed that for strong magnetic fields, if the star 
surface temperature is below a critical value (which is a 
function of the magnetic field strength), the atmosphere 
can undergo a phase transition into a condensed state 
(see also [Lai 2001). For B > lO^'' G, this may occur 
even for temperatures as high as 10® K. This raises the 
possibility that the thermal radiation is emitted directly 
from the metal surface of the NS. 

The thermal emission from condensed F e surface of 
magne tic NSs was previously studied by [B rinkman^ 
10m (see also l!tohlll975tlLenzen fc TriimDeri.1978.) and 
shown to produce a rough blackbody with reduced emis- 
sivity and a spectral feature at the electron plasma en- 
ergy. For the temperatures and magnetic fields (T > 
10^ K and B — 10^^~^^G, appropriate for accreting X- 
ray pulsars) considered by Brinkmann, the Fe surface is 
not expected to be in the condensed state. However, 
at lower temperatures appropriate for dim, isolated NSs, 
or for higher B appropriate for magnetars, condensation 
remains a possibility (see Lai 2OOI4). 

In this paper, motivated by recent observations of dim 
isolated NSs, we calculate the emissivity of condensed Fe 
or H surface of magnetic NSs in the regime where we 
expect condensation mi ght be possible. Our study goes 
beyon d previous work IjBrinkmannl 1198(1 iT^irolla et al.l 
I2004D in that we calculate both the spectrum and po- 
larization of the emission, and provide a more accurate 
treatment of the dissipative effect and transmitted radi- 
ation. In previous works, the ions have been treated as 
fixed; while the exact dielectric tensor of the condensed 
matter is currently unknown, we also consider the alter- 
nate limit of free ions (see §2.2). 

Regardless how the effect of ions in the dielectric ten- 
sor is treated, we find appreciable difference between our 
result and that of TuroUa et al. We traced the difference 
to their neglect of the ion effect, and their "one- mode" 
treatment of the transmitted radiation in the low-energy 
regime (see §4.1). Some of our preliminary results were 
reported in I Arras fc Lail 1)19991) . 

This paper is organized as follows. Section |21 sum- 
marizes the basic properties of the condensed matter in 
strong magnetic fields. The method for calculating the 



emission from the surface is outlined in ^ and numeri- 
cal results presented in ^ We discuss the implications 
of our results to observations of dim isolated NSs and 
magnetars in ^ 

2. CONDENSED SURFACE OF MAGNETIC NEUTRON 
STARS 

2.1. Condition for Condensation 

It is well known that strong magnetic fields can quali- 
tatively change the properties of atoms, molecules and 



condensed matter. For B ^ Bq 



Z 



2^3, 



2/h^ 



2.35 Z^ X 10^ G (where Z is the nuclear charge num 
ber), the electrons in an atom are confined to the ground 
Landau level, and the atom is elongated, with greatly 
enhanced binding energy. Covalent bonding between 
atoms leads to linear molecular chains, and interactions 
between molecular chains can lead to t he format ion of 
three-dimensional condensed matter fsee lLaill2001l for a 
recent review). 

For H, the pha se diagram und e r diff erent conditions 
has been studied. iLai fc Saloeteil l)1997|) showed that in 
strong magnetic fields, there exists a critical tempera- 
ture Tcrit below which a phase transition from gaseous 
to condensed state occurs, with kTcrit about 10% of the 
cohesive energy of the condensed hydrogen. Thus, Tcrit '^ 
8 X 10^, 5 X 10^ 10® K for B = lO^^, lO", 5 x 10^* G 
l)Laill200l() . An analogous "plasma phase transition" was 
also obtained in an alternative thermod ynamic mode l 
for magnetized hydrogen plasma ( Potckhi n et al.ll999al). 
While this model is more restricted than Lai & SalpcteJ 
L1B21) in that it does not include long H„ chains, it treats 
more rigorously atomic motion across the strong B field 
and Coulomb plasma nonideality. In the Potekhin et al. 
model, th e density of phase sep aration is roughly the 
same as in ILai fc SalneteH ^9T) (see eq. P below), but 
the critical temperature is several times higher. Thus 
there is probably a factor of a few uncertainty in Tcrit. 
However, there is no question that for T < Tcrit/2, the H 
surface of the NS is in the form of the condensed metallic 
state, with negligible vapor above it. 

For heavy elements such as Fe, no such systematic 
characterization of the phase diagram has been per- 
formed. Calculations so far have shown that at 10^^ — 
10^"^ G, a linear chain is unbound re l ative to individ- 
ual a toms for Z > 6 Ij.TonesI 119861: INeuhauser et al 
1987) - contrary to earlier expectations (Flowers et al 
1977(1 .® Therefore chain-chain interactions play a cru- 



cial role in determining whether 3D zero-pressure con- 
densed matter is bound or not. Numerical results of 
Ijonegj (|1986j) . together with approximate scaling rela- 
tions suggest an upper limit of the cohesive energy (for 

Z > 10) Qs < Z^/'^B%^ eV, where B12 = B/{10^^ G). 
Thus for Fe, the critical temperature for phase transition 

Tcrit < O.lQs/k < lO^-^Sif K (Lai 2001!) . 

The zero-pressure density of the condensed matter can 
be estimated as 



Ps ~ 560 r; ^ Z-^/^ B%^ g cm 



-3 



(1) 

where A is the mass number of the ion {A w 1.007 for 
H, ^ w 55.9 for Fe), and 7/ = 1 corresponds to the uni- 

6 For sufficiently large B, when B > W^{Z/2&j^ G, we expect 
the linear chain to be bound in a manner similar to the H chain 
lLaai2001]) . 
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form electron gas mod el in the Wigner-Seitz approxima- 
tion IjKadomtsevll 1 97flD . Other effects (e.g., Coulomb ex- 
change interaction, nonuniformity of the electron gas) 
can reduce the density by up to a fac t or of ~ 2, and 
thus 77 may be as small as 0.5 l|Lail 12001"; see also 
iPote khin & Chabrieill20M) . The condensate will be in 
the liquid state when the Coulomb coupling parameter 
r = {Zef/{a^kT) = Q.227 Z'^{pi/AY'^ /Tq < r„. Here, 
Qi is the ion sphere radius ((47ra^/3)~^ = n^, where rii is 
the number density of ions), pi = Ps/{^ g cm~^), Tq = 
r/(10^ K), and r^ is the characteristic value of F at 
which the Coulomb crystal melts. In the one-component 
plasma model (i.e., classical ions on the background of 
the uniform degenerate electron gas), r,„ — 175, but the 
electron gas nonuniformity (i.e., electron screening) in- 
troduces a dependence of T m on p and Z, typically withi n 
the range T^ -- 160-190 (jPotekhin fc Chabried l200(l . 

From equation (jTJ we obtain F ~ 1.876 t]^^^ Z^/^ B^^^ /Tq 
at the condensed surface. Therefore, the surface will be 

solid when T < 7 x 104ryi/3(175/F,„)Si{^ K for H (where 



,2/5 



Si4 = B/10'^ G) and T < 4 X 10'^?7^/^(175/r„0Bi^" K 
for Fe. Therefore, if condensation occurs (T < Tcrit), 
we expect the Fe condensate to be solid. Note that we 
use the simple melting criterion above for the condensed 
phase only. It cannot be used for non-condensed iron 
at T < 10'^ K (e.g., when T is only slightly above Tent) 
because in this case the state of matter is affected by 
partial ionization. 



2.2. Dielectric Tensor of Condensed Matter 

The emissivity of the condensed NS surface will depend 
on its (complex) dielectric tensor (see 33 • As a first ap- 
proximation, we consider the free electron gas model fo r 
the condensed matter (e.g., Ashcro ft fc Merm in '1976^. 
In the coordinate system with magnetic field B along the 
z-axi s, the dielectric tensor takes the form (cf. iGinzburel 

ri97nli 7 



[-] 



z=B 




(2) 



where 



e±5~l- 



r/~ 1 



(l±uy')(lTu, 



V2,^^^(*.) 



1 



I 7 



«■ 



(3a) 
(3b) 



In eq. ij^j, the dimensionless quantities Ue = {Ebc/E)'^, 
Ui = {Ebi/EY, Ve = (Epe/E)^ are used, where E = huj 
is the photon energy, Ebs, Ebi are the electron and ion 
cyclotron energies, and Ep^ is the electron plasma energy. 



■^ See also lLai fc Hcl2nn3al Note that eq. (13) of lLai fc Hcl2nn3al 
is incorrect: 7^^ should simply be 7ei(l + Zme/Arrip). We neglect 
the factor 1 -|- Zm^/Arrip in eq. I3ai since it provides a negligible 
correction relative to the uncertainty in the coUisional damping 
(see §2.3). 



These energies take the values: 
EBe = — = 1158^14 keV, 



Ep 



En 



nieC 
hZeB 



0.635 Si4 -r keV, 
mic \ A , 

. ^o o X 1/2 



Z 



(4a) 
(4b) 



1/2 



1/2 



0.0288 ( - ) pl' keV 



= 10.8 7/1/2^1/^ B^f keV, 



(4c) 



where n^ is the electron number density, and nii is the ion 
mass. The coUisional damping is calculated for motions 
transverse and longitudinal with respect to the magnetic 

(tr) (I) 

field. The dimensionless damping rates 7^^ and 7^/ are 
obtained from the coUisional damping rates v^,/" and 1^1,^ 



(see 3131) through j^p = Hiyif^/E and 7^'^ = hiy^^/E. 

Equations ^ give the elements of the dielectric ten- 
sor for a cold, magnetized plasma. While the expres- 
sions were derived classically, the quantum calculation, 
incorporating the quantized nature of electron motion 
trans verse to the magne tic field, vields identical results 
(e.g. JCanuto fc Venturai.l972:.Pa^v et alJll980|) . More 
significantly, expressions (O assume that the electrons 
and ions are subject to the pairwise Coulomb attraction, 
the interaction with the stationary magnetic field, and 
the periodic force from the propagating electromagnetic 
wave. At high densities, however, other interactions can 
also be important. For instance, the ions are strongly 
coupled to each other when the Coulomb parameter F 
is large. It is this coupling that leads to the liquid- 
solid phase transition mentioned in §2.1. One might 
suggest that in the solid phase the ion motion should 
be frozen (by se tting the ion mass ru j — 00), as implic- 
itly adopted bv iTuroUa et all l|2004|) . This is not ex- 
actly true. It is known that optical modes of a crys- 
tal lattice (at B — 0) can be described by polarizabil- 
ity of the form given by equation ^ with an additional 
term in the deno minator which specifies the binding of 
the ions (see, e.g.. IZimanlll979D . Acco rding to th e har- 
monic model of the Coulomb crystal IjChabried 119931 . 
the characteristic ion oscillation frequency (the Debye 
frequency of acoustic phonons) is ujo ~ OAEpi/h, where 

Epi — 6.75 X 10^* {Z/A) p^' keV is the ion plasma en- 
ergy. The magnetic field appreciably affects the motion 
of the ions in the Coulomb cry stal if Hlub/Ebi ^ 1 (or 
Epi < Ebi, see lUsov et al.lll980|l . From equations Q) we 
find hwo/EBr « 1.6r;i/2^i/22'-0-3B-0-4^ which shows 
that the magnetic forces on the ions are not completely 
negligible compared to the Coulomb lattice forces. 

Needless to say, our current understanding of the con- 
densed matter in strong magnetic fields is crude, and 
equations |(2J) are only a first approximation to the true 
dielectric tensor of the magnetized medium. In our cal- 
culations below, in addition to the case of of quasi-free 
ions described by equations l(2Jl, we shall also consider the 
case where the motion of the ions is neglected (formally 
obtained by setting mi — 00). It is reasonable to expect 
that in reality the surface radiation spectra lie between 
the results obtained for these two limiting cases. Never- 
theless, future work is needed to evaluate the reliability 
of our results at low frequencies. 
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Fig. 1. — Transverse and longitudinal damping rates 7l!/^^ ' and 

hu^^ as a function of magnetic field strength B = 10^^i?i2 for 
condensed Fe surface at T = 10^ K. The density is calculated 



using eq. 0, with rj = 1. 



2.3. Collisional Damping Rate in the Condensed Matter 

(I tr) 

For the collisional damping rates 7^/ , different ap- 
proximations can be used in different ranges of frequency 
CO and density p. For E ^ Epe = hujpe (where ujpe 
is the electron plasma frequency), the electron- ion col- 
lisions can be considered as independent, and 7^/ ^ are 
determined by the effective rates of free-free transitions 
of a single electron-ion pair. However, this approxima- 
tion fails at -E < Ep^, where collective effects become 
important. Moreover, the electron degeneracy should 
be taken into account in the condensed surface. In 
general, the complex dielectric tensor e for arbitrary w 
can be obtai ned from kinetic theory, at least in princi- 
ple (e.g., see lGinzburelll970|) . Since such a general ex- 
pression of e is unknown at present, we shall approxi- 
mate 7gj' ^' in the E < Epe regime using the result of 
Potekhin (1999), who obtained the zero-frequency con- 
ductivity tensor for degenerate Coulomb plasmas (liq- 
uid and solid) in arbitrary magnetic fields. Specifically, 

we set i/gj' = 1/t||, i/^/ = 1/t±, where r|| and t± 
are the effective c ollisio n times given by eqs. (28) and 
(39) of iPotekhml l)1999j) . respectively. Figure Q] shows 



fi-z/g/ and fifgj as a function of magnetic field strength 
for condensed Fe surface at T = 10^ K, over the range 
B = IO12 - 10" G. 

The calculations of i^^/ and z/^/ adopted in our paper 
neglect the influence of the magnetic field on the motion 
of the ions. Therefore, these calculations apply only in 
the Mi ^ limit (this corresponds to the "fixed" ion 
limit of §2.2), or in the regime E ^ Est- We note, 
however, that the emissivity at E < Ebi does not depend 
sensitively on the damping rates (see §4; in particular. 
Fig. 2 shows that the emissivity at such low energies is 
almost the same with/ without damping). Thus, unless 



orders of magnitude larger than our adopted values, our 
emissivity results will not be affected by this uncertainty 
(Indeed, as discussed in §2.2, the main uncertainty at 
such low energies lies in whether to treat the ions as 
"free" or "fixed"). 

3. EMISSION FROM CONDENSED MATTER: METHOD 

In this paper we consider the regime where a clear 
phase separation occurs at the NS surface (i.e., for T 
at least a few times lower than Tcrit), so that the vapor 
(gas) above the condensed surface has negligible density 
and optical depth. In this case the radiation emerges 
directly from the condensed matter. 

3.1. Kirchhoff's Law For A Macroscopic Object 
A macroscopic body at temperature T produces an in- 
trinsic thermal emission, with specific intensity I^ ■ To 
calculate the intensity, consider the body placed inside 
a blackbody cavity also at temperature T, i.e., the body 
is in thermodynamical equilibrium with the surrounding 
radiation field, whose intensity is given by the Planck 
function B^{T). Imagine a ray of the cavity radiation 
impinging on a surface element dA of the body. The ra- 
diation field is unpolarized, and the electric field of the 
incoming ray can be written in terms of two indepen- 
dent polarization states: E-[ — Ae^ and Ej — ^63 , 

where A = yjByjl^ and &\ and 62 are the polarization 
eigenvectors of the incident wave. The ray is, in general, 
partially refiected, each incoming polarization giving rise 
to a reflected field: 



eI'") = A 



E, 



(r) 



= A 



(r) 

(r) 
^2161 



(r) 
(r) 



(5a) 
(5b) 



.(r) 



.(r) 



where E]^ ' and E2 are the reflected electric fields due 

to incoming fields e![ and E2 , respectively. Thus, the 
intensity of radiation in the reflected fleld with polariza- 



ions e^ and 62 is: 






4f^-^(|mr + hi|>. 


— l^RlBiy, 


(6a) 


4? = i(ki2r + |r22r)i?. 


^ 2^2^.. 


(6b) 



The energy in the incoming wave for frequency band 
1/ ^ I' + diy during time dt is B^ dAdH,^^' dvdt, where 
dfi^*-' is the solid angle element around the direc- 
tion of the incoming ray. Similarly, the energy con- 
tained in the reflected wave (along each polarization) 
is {l/2)Ri^2B^dAdn^''Ui^dt, with df^W = df^W. To 
insure that the cavity radiation fleld remains an unpo- 
larized blackbody, the intensities of radiation emitted by 
the body (in the same direction as the reflected wave) 

(r) (r) 

with polarizations e\ and ej must be: 



r(e) _ 



Ar) 



r(e) 



r(e) 



(7a) 
(7b) 



the true values of v. 



(0 .,(*'■) 



Since I^i and 1^,2 are intrinsic properties of the body, 
at such low energies are many these equations should also apply even when the body 
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is not in thermodynamical equilibrium with a blackbody 
radiation field. Thus, a body at temperature T has emis- 
sion intensity 

li^'^ = {l-R)B,{T) = JB,{T) (8) 

where R = (1/2) (i?i + R2) is the reflectivity, and J — 
1 — i? is the dimensionless emissivity. The degree of hnear 
polarization of the emitted radiation is 



r(e) _ r(e) 



/ 



(e) 



(e) 



1 R2 — Ri 

2 1- R ' 



(9) 



In the xyz coordinate system, the rotated dielectric 
tensor takes the form: 

ecos^ Ob+t] sm^ 9b igcosOB {e — ri)sm 6b cos 9b 
e|=l —ig cos 9b e — igsin^s 

, (e — 77) sin 9b cos 9b ig sin 9b e sin^ 9b + V cos^ 9bJ 

(13) 
For eq. (|12|l to have a non-trivial solution, the deter- 
minant of the matrix e + {n- '^ ( k^- k^- — ij must be 
equal to zero. This gives an equation involving powers of 



3.2. Calculation of Reflectivity 

To calculate the reflectivity i?, we set up a coor- 
dinate system as follows: the surface lies in the xy 
plane with the z-axis along the surface normal. The 
external magnetic field B lies in the xz plane, with 
B X z = sin^sy, where 9b is the angle between B 

or 



(*) 



, sin 9, , and cos 6':- . Substituting eq. Ullfl into this 
equation, and squaring both sides yields a fourth-order 
polynomial in {n, )^, which allows for the determination 
of the indices of refraction (see Appendix IXI for more de- 
tails). Having determined n- , eq. ((T^ can be used to 



lit) Jt) 



(0 



AD • z = 0, AB ■ z = 0, AE X z = 0, AH x z = 0, 

(14) 
where, e.g., AE = EW -^E^ - E^, D^*) =e-E(*), and 



calculate e^- , while eq. (fTT|l gives 9^ . Once 0^- , e 
and z. Consider a ray (of certain polarization, e^*^ or and n^ are known, rn, ri2, in and ti2 can be obtained 

e^'^) impinging on the surface, with incident angle 6'(') using the standard electromagnetic boundary conditions: 

and azimuthal angle cp, so that the unit wave vector 
kW = (-sin6'Wcos(/?, -sin6'Wsinv?, -cos^W). The 
transmitted (refracted) and reflected rays lie in the same 
plane as the incident ray. Our goal is to calculate the 
field associated with the reflected ray. 

Outside the condensed medium {z > 0), the dielec- 
tric tensor and permeability tensor are determined by 
the vacuum polarization effect with e = al + qBB and 
fj, ^ = al + Tn BB, where o,q.m are dimensionless func- 
tions of B fsee lHo fc Laill200^ and references therein). 
Since a '- 1 and |g|, |m| < 1 for B < 5 x 10^'' G, the 
vacuum polarization effect is negligible. In our calcula- 



HW=B(*)=^ ni^'ink^xei' 



,(*) 



4h,2ki'^ X e^i^ 



tion (Appendix A), we choose e]^' (and e]^') to be along 



J«) 



,W^ 



the incident plane and 63 (and 63 '' ) perpendicular to 
it. 

Consider an incident ray with E^*' = E]^'^ = Ae]^. 
The E-field of the reflected ray takes the form given by 
eq. ()5a(l . while the transmitted wave has the form: 



EW=Ef'=^ tnei*^+ii2e^' 



,(*) 



(10) 



The eigenvectors of the transmitted wave, e\ ' and 63 



(t) 



.(*) 



(«) 



3(t) 



depend on the refraction angles, 9\ and ^3 j respec- 
tively; note that in general, these angles are complex 
and different from each other. The refraction angle 9, , 

the mode eigenvector e^- and the corresponding index 

of refraction n, (j = 1, 2) satisfy Snell's law 





sin6'W =n^.*^sin6' 




and the mode equation* 




e-f(nf)^(kfkf -l)' 


• e(*) = 0, 


where I is the unit tensor, and kj*^ 


{-sm9f 


cos 1^, — sm d)j sm (^, — cos 


i9f) is the 



(11) 



(12) 



unit 



(15) 
neglecting the vacuum polarization effect (/i, ~ I). Note 
that eqs. 1)14(1 are not all independent. Only AE x z = 
and AB x z = are used in our calculation. 

A similar procedure applies in the case when the inci- 
dent wave is E^'^ = E2 = •^^2^ yielding the reflection 
coefficients r2i and r22 (together with ^21, ^22)- 

4. EMISSION FROM CONDENSED SURFACE: RESULTS 

In this section, we present results of surface emission 
for three illustrative cases: Fe surface at B = 10^^ G and 
10^3 G, and H surface at W'^ G. As discussed in EH 
the condensation temperature for these cases is around 
10^ K. Note that the dimensionless emissivity J = 1 — 
R [see eq. ©] depends only weakly on T through the 
collisional damping rate Ci'Z.'Zy For concreteness, we set 
T = 10® K in all our calculations. Figures |301 show the 
emissivity J as a function of photon energy E in the three 
cases; the B fleld is assumed to be normal to the surface 
{9b = 0). In all three cases, the emissivity is reduced 
from blackbody at low energies, while approaching unity 
for E > a few x Ep^. In the case of Fe, there are features 
associated with the ion cyclotron energy Ebi and the 
electron plasma energy Ep^. For H, the electron plasma 
energy is too high to be of interest for observation, but 
the feature around the ion cyclotron energy is evident. 

The spectral feature in the emissivity J near Ebi can 
be understood by considering the special case of 6'(') — 
(normal incidence). In this case the reflectivity takes the 
analytic form: 



wavevector of the transmitted waves. 



* The vacuum polarization effect is neglected in eq. 1121 . 
which is justified because the density of the condensed medium 
is much larger than the vacuum resonance den sity, py ~ 
0.96(A/Z)Bl^(E/keY)^ g cm-3 (see lLai &: H^l200a . 
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where ni and 712 are the indices of refraction of the two 



modes in the medium, and are given by nj 
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Fig. 2. — Dimensionless emissivity J = 1 — R as a, function of 
photon energy E (keV) for the case of condensed Fe surface, at 
B = 10^^ G. The B field is normal to the surface. The different 
curves correspond to different angles 9'*' between incident photon 
direction and surface normal. The short-dashed-dottcd line (la- 
beled "no damping") shows the result when the coUisional damp- 
ing is set to zero in the plasma dielectric tensor. The other light 
lines (labeled "no ion") show results when ion motion is neglected 
for two values of 6'*' (i.e., by setting the ion mass to oo; see §2.2). 
The three vertical lines denote the ion cyclotron ene rgy Egi, the 
electron plasma energy Epe [see eq. jl}] and Ec [eq. 1191 ]. 
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Fig. 3.— Same as Fig. 2, except for B = lO" G. 
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Although 7gj'^-^ can be greater than unity (see Fig.^, the 
imaginary part of nl 2 can be neglected for a qualitative 



0.8 



e(i 

. e<i 



= 
= 71/6 

= tt/4 

= 71/3 

- (no ion) 

- tt/3 (no ion) 



0.4 



0.2 







1 


1 


0.001 


0.01 


0.1 








E (keV) 



r 



Fig. 4. — Same as Fig. 2, except for H surface at 10^ G. 



understanding of the spectral features, since Ve/ue <IC 1 
[see eq. (0J]. Then eq. l(T7|) becomes 



'1,2 



It 



1/2,1 l/2\ 



(18) 



For E < Esi {ui > 1), both ni and 712 are real and 
differ from unity, leading to J < 1. For E > Egi, ni is 

imaginary until {ve/uj )(1 — u,-J )^^ < 1, which occurs 
at 

Ec^Ee^ + ElJEBe- (19) 



Thus, for Eb, < E < E, 



c, n\ 



increases from — cxd to 



(implying no mode propagation in the medium), giving 
rise to the broad depression in J (with J — > 0.5 as the 
energy nears Ec)- 

We can similarly understand the feature near the elec- 
tron plasma energy. This feature appears only for 9^"^' ^ 
0. For energies around Epe, Ue ^ 1, Wi <c; 1, and we have 

1/2 

e ~ 1 + Ve/Uf. and g ~ —Ve/ue ■ Substituting these val- 
ues into H12|l and neglecting terms to order fe/ue and 
higher, we find 

6*^, nl^l (20) 



i^-Ve) 



For E > Epe, both ni and n2 are real, while for E < Epe, 
n\ < 0. The reflectivity no longer takes the simple an- 
alytic form of (|16|l . However, the basic behavior of the 
reflectivity is largely the same: for one mode with imag- 
inary n and the other with n w 1, the emissivity J at- 
tains a local minimum (~ 0.5 in the absence of collisional 
damping; see Fig. [21). 

When calculating the emissivity, it is clear that the in- 
clusion of the ion terms in eqs. (PJ for the elements of the 
dielectric tensor can qualitiatively change the emission 
spectrum at low energies (see Figs. |21[Z|)- As discussed 
in §2.2, complete neglect of ion effects is not justified; 
while the exact dielectric tensor is currently unknown, 
the true spectra should lie between the two limiting cases 
we present here. Without the ion terms, the broad fea- 
ture around Ebi is replaced by a stronger depression of 
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Fig. 5. — Dimensionless emissivity J = 1 — R as a, function of 
photon energy E for the case of condensed Fe surface at _B = 10^'' 
G. The incident angles are fixed at 9'*' = 7r/4 and <p = 7r/4. The 
different curves correspond to different magnetic field inclination 
angles [Bg is the angle between B and the surface normal). As 
in Fig. 2, the light lines (labeled "no ion") show results when ion 
motion is neglected in the plasma dielectric tensor. 
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Fig. 7. — Spectral flux as a function of photon energy E for the 
cases of condensed Fe (B = lO^^ , lo" G) and H (B = lO" G) 
surfaces, all at temperature T = 10^ K. The light lines (labeled 
"no ion") show the flux for Fe and H surfaces when ion motion is 
neglected. The solid line shows the blackbody spectrum at lO" K. 
For all of the curves, the magnetic inclination angle Bb = 0. 
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Fig. 6. — Same as Fig. 5, except that the geometry is fixed at 
0(0 = 7r/4 and dg = 7r/4, and the different curves correspond to 
different (/3 (t he an gle of the plane of incidence with respect to the 
xz plane; see i|3.2l . 



J at low energies, up to i? ^ Ec- At high energies, the 
ion effect is unimportant. 

Figures |S1 and give some examples of our numerical 
results for the cases when the magnetic field is not per- 
pendicular to the surface (i.e., 6b t^ 0). In these cases the 
emissivity J is no longer symmetric with respect to the 
surface normal, but depends on 6^'^' and the azimuthal 
angle (/s. Although the geometry is more complicated, 
the basic features of the emissivity are similar to those 
depicted in Figs.|21|3| 

Figure |7| depicts specific flux at the NS surface, 
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Fig. 8. — Degree of linear polarization P [see eq. J5J] as a func- 
tion of photon energy E for the case of condensed Fe surface, with 
B = 10^"^ G. The B field is normal to the surface, and the different 
curves correspond to different angles 0'^' between incident photon 
direction and surface normal. The net linear polarization is peaked 
around Egi and Epe- Positive P corresponds to polarization paral- 
lel to the k-B plane, while negative P corresponds to polarization 
perpendicular the k-B plane. Note that P changes sign around 
Ec- 



F^ = /„^''d(^/J'^^d6'(*)cos6'«sin6iWj(6iW,(p)B^(T), as 
a function of photon energy for the three cases illustrated 
in Figs.|2JS For the Fe surface, there is a reduced emis- 
sion (by a factor of 2 or so) around Esi ^ E < Ec com- 
pared to the blackbody at the same temperature. For the 
H surface at B = 10^"' G, the flux is close to blackbody 
at all energies except for a broad feature around Ebi- 
The radiation from the condensed surface is polarized. 
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Fig. 9. — Same as Fig. 8, except for the H surface at B = lO^** G. 
There is a slight net hnear polarization perpendicular to the k-B 
plane (P ~ —5%), except around Esi where the polarization peaks 
parallel to the k-B plane. 
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Fig. 10. — Angle-averaged intensity (1 — i?) as defined in §4.1 for 
_B = 5 X 10^3 G, Te = 1.0, Qb = 0.7 x 7r/2. The solid line shows 
our result including the ion eff'ect, while the dashed-line shows the 
results when the ion motion is neglected. For comparison, the 
dotted line shows data from Fig. 5 of Turolla et al. (2004). 



Figures |H1 and |3| show the degree of hnear polarization 
as a function of energy for the cases iUustrated in Fig- 
ures Eland |3|(i.e., Fe at 10^3 Q and H at 10^'' G). The 
degree of hnear polarization P increases with angle of in- 
cidence, and is clearly peaked around Ebi and Ep^ . For 
the Fe surface, at energies below Ec^ the polarization 
vector is parallel to the k-B plane. Above Ec, the sign 
of P changes, and the radiation is polarized perpendicu- 
lar to the k-B plane. For H, there is a slight net linear 
polarization perpendicular to the k-B plane, except near 
EBi^ where the polarization peaks with P > 0. These 
polarization properties of condensed surface emission are 
qua litatively different from those of atmosphere emission 
('see lLai fc Hdl20r)3U and references therein). 

4.1. Comparison with Previous Work 

Recently, Turolla et al.(2004) performed a detailed cal- 
culation of the emissivity of a solid Fe surface. Our re- 
sults differ significantly from theirs in several respects. 
In particular, Turolla et al. found that collisional damp- 
ing in the condensed matter leads to a sharp cut-off in 
the emission at low photon energies, especially when the 
magnetic field is inclined with respect to the surface nor- 
mal. For comparison, in Fig. 1101 we show the angle- 
averaged emissivity, (l — R) = F^/[TrB,^{T)], for a specific 
case with B = 5x lO^^ q y^ = ^ q and 9b = 0.7 x tt/2; 
this should be directly compared with Fig. 5 in Turolla et 
al. Their results show no emission below ~ 0.1 keV, and 
they find that this "cutoff" feature becomes more pro- 
nounced as the magnetic field inclination angle increases 
and the field strength decreases. Our calculations clearly 
do not show this behavior (see the sohd line of Fig. ITIHl . 

These discrepancies stem from at least two differences 
in the reflectivity calculation: (1) Turolla et al. neglected 
the effect of ion motion in their expression for the plasma 
dielectric tensor (see the end of §2.2). This strongly af- 
fects the emissivity a.t E < Ebi (see also Figs. |2[7I). 
(2) Even when the ion motion is neglected (by setting 



nii — cxd), our result (see the dashed-line in Fig. 10) 
does not reveal any low-energy cutoff. It is most likely 
that this difference arises from the "one-mode" descrip- 
tion for the transmitted radiation adopted by Turolla et 
al.: when the real part of the index of refraction of a 
mode is less than zero or the imaginary part of the in- 
dex of refraction exceeds a threshold value, this mode 
is neglected by Turolla et al. in the transmitted wave. 
Such treatment is incorrect and can lead to significant 
errors in the reflectivity calculation. The inclusion of 
collisional damping gives rise to complex values for the 
index of refraction, which lead to transmitted waves with 
a propagating (oscillatory) part multiplied by a decaying 
amplitude (see Appendix B). While the damping factor 
for such waves can be large if the index of refraction has 
a large imaginary part, the propagating piece allows en- 
ergy to be carried across the vacuum-surface boundary, 
therefore these waves cannot be ignored in the reflectivity 
calculation.^ 

5. DISGUSSION 

As discussed in §1, many isolated NSs detected in ther- 
mal emission display no spectral features and are well fit 
by a blackbody spectrum. The most thoroughly studied 
object of this type is RX J1856.5— 3754, which is well fit 
in the X-ray by a blackbody spectrum at A: Too = 63.5 eV, 
with emission radius Roo = 4.4 (d/120 pc) km (where d 
is the distance). This X-ray b lackbody und erpredicts 
the optical fiux by a factor of 7. iPavlov fc Zav lin ( 200v|) 
review several models involving a non-uniform tempera- 
ture on the surface of the NS, in which the X-ray photons 
are emitted by a hot spot. By varying the temperature 
distribution and assuming blackbody emission from each 

^ Af ter submitting our paper to ApJ, a preprint by 
IPerez-A zorin et al. 1 2004} appeared on the archive, reproducing 
the calculations described here. They arrive at conclusions similar 
to the ones discussed above. 



Radiation from Condensed Surface of Magnetic Neutron Stars 



surface element, a reasonable fit to the X-ray and opti- 
cal data can be achie ved (see also lBraie fc R.omanill2002l: 
iTriimper et alJl2003() . Nevertheless, the nearly perfect 
X-ray blackbody spectrum of RX J 1856. 5— 3754 is sur- 
prising. 

If the NS surface is indeed in the condensed form (see 
§2), the emissivity will be determined by the properties 
of the condensed matter. Our calculations (§3 and §4) 
show that the emission spectrum resembles that of a di- 
luted blackbody, with the reduction factor in the range 
of J = 0.4 — 1 depending on the photon energy (see Figs. 
[2Hnil- This would increase the inferred emission radius 
by a factor of J~^^^ . The weak "absorption" features in 
the emission spectrum are associated ion cyclotron res- 
onance and electron plasma frequency in the condensed 
medium. We note that the emissivity and spectrum pre- 
sented in this paper correspond to a local patch of the 
NS. When the emission from different surface elements 
are combined to form a synthetic spectrum, these ab- 
sorption features are expected to be smoothed out fur- 
ther because of the magnetic field variation across the 
NS surface. 

In our calculations, we have assumed a perfectly 
smooth surface. This is valid if the condensed matter 
is in a liquid state, as is likely to be the case for H 
condensate (see §2.1). For Fe, the condensed surface 
is most likely a solid and we may expect a rough sur- 
face. Although it is not possible to predict the scale 
and shape of the surface irregularities, their maximum 
possible height /imax can be estimated from the require- 
ment that the stress nonuniformity '^ pgh is small com- 
pared to the shear stress ii Q^- With the shear modu- 
lus ii ~ 0.1 n^ {ZeY /a^ (e.g., JOgata fc IchimanJ (|199nr >) 
and the maximum strain angle 9s = 1O~^0_3, we find 

/imax -^ 2 X 10-5 0-3 Z'^ A'^l'^ p^'^ g^l cm (where g = 
VS^^g\4, cm s"^ is the surface gravity). For the con- 
densed Fe surface at the density given by eq. (Q), we 

have /imax ^ 4 X 10^^ d-z B^2 cm (for a neutron star 
with i? = 10 km and M = IAMq). Clearly, the scale of 
the surface roughness can easily be much larger than the 
photon wavelength ('^ lOA). As illustrated in Fig. IllL 
the surface may be much less reflective than the results 
shown in §4, and the emission will be closer to the black- 
body spectrum. 

The emission from a condensed NS surface is distinct 
from atmospheric emission in several aspects: (i) Atmo- 
spheric emission generally possesses a hard spectral tail 
(although this tail is somewhat suppressed by the QED 
effect for S > 10^^ G; see §1), whereas the condensed sur- 
face emission does not; (ii) The spectrum of a cool NS 
atmosphere can have both cyclotron and atomic absorp- 



tion features (again, they are reduced for B > 10^"^ G); 
the broad (cyclotron and plasma) features of condensed 
surface emission persist even in the magnetar field regime 
(if they are not smoothed out by variations of surface B 
fields or by the rough surface effect); (iii) The polariza- 
tion signature of condensed matter emission is qualita- 
tively different from that of atmospheric emission. All 
these aspects can serve as diagnostics for the physical 
condition of the emission region. 

At the surface temperature of AXPs and SGRs, T ~ 
5 X 10^ K, H is unlikely to be condensed, but Fe con- 
densation is possible. The dim, isolated NSs have lower 




Fig. 11. — Effect of surface roughness on the reflectivity. The 
surface roughness is characterized by the vertical scale h and hor- 
izontal scale I, both much greater than the photon wavelength. 
For the idealized "triangular" surface, a normal incident ray goes 
through at least two reflections if 9 = tan~^{i//i) < 60°, at least 
three reflections \{ 9 < 36°, at least four reflections if S < 180° /7, 
etc. Thus net reflectivity of the rough surface is ^ 1 if /i > a fewxi, 
and the emission spectrum will be close to blackbody. 



temperatures (T < 10^ K), and if they possess magnetar- 
like fields, condensation is likely. In particular, the black- 
body X-ray spectra of RX J1856. 5-3754 [kT - 64eV) 
and RX J0420-5022 {kT ~ 45eV; see iHaberl et all 
[2004b) could arise from condensed surface emission (e.g., 
non-smooth Fe surface at i? > 10^^ G), although to ac- 
count for the optical data, nonuniform surface tempera- 
tures are still needed. 
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APPENDIX 



A. REFLECTIVITY CALCULATION 
Here we fill in some of the details for the reflectivity calculation described in ti3.2l 
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In the coordinate system xyz defined in M3.2I the explicit expression for eq. (|12|l is: 
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where Uj (with _7 = 1, 2) is the index of refraction in the medium, and 9a is the formal complex angle of propagation 

calculated using Snell's law (see Appendix B for a discussion of the interpretation of complex 6*; ). Taking the 
determinant of eq. I|A1|I yields: 



a^n^ + 02^^- + cos6'j- sin6'''*'' {a\nj + aa'^j) + oo = 0, 
where we have used Snell's law and the following definitions: 
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The cos 9) term is moved to the right-hand side, and the entire equation is then squared. Using the identity 
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COS 6*^- = 1 — sin 9j and Snell's law gives a polynomial equation in 
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(A4) 

The polynomial equation (|A4p has eight roots for rij which we found numericall y us ing Laguerre's method (Press et 
al, 1996). Only two of the roots are physical and satisfy the original equation (|A2I) . In practice, it was found that 
for certain combinations of the parameters E, 0(*\ 9b, (p, a spurious root satisfies eqn. ljA2p to the specified degree of 
accuracy, resulting in an unphy sica l resul t fo r the reflectivity. It is often the case that such roots can be discounted 
physically using the conditions (|B3(I and ljB4|l (see Appendix B). Once the indices of refraction tt-i, n2 are known, the 

normal mode polarization vectors can be determined. Solving eqn. (|12(l for the ratios Ex /Ey and Ez / Ey gives 
the expressions: 
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With the propagation modes in the plasma determined, the latter two equations of ifTHl give: 
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sm ip 



/coseiW 
cos 9^''^ cos If 

— COS(p 

\ — sin (fi 



COS if Ci C2 \ 

simp — Ci — C2 

— cos 6*^'^ sin ip C5 Cq 

cos 6*^ cos ^ C7 Cs ) 




— cosfy^'^ sm(y9 
cos 0^'^ cos Lp 

C0Sly9 
SVi\p> 

— cosip 

— SVHip 

— cos0^') siniy9 

cos 0^'^ cos Lp 

.(■0 - 



(A6) 



for the incoming polarization modes e^ = (— cos^'^'^ COS93, — cos0^'^ siniy9,sin6''^'^) and 62 — (sin(p, — COS93, 0). In- 
verting the coefficient matrix of eqn. (|A6|) and performing extensive algebra yields the following expressions for the 
reflected field amplitudes: 



r\\ 



ri2 



r2i 



r22 



4A cos e'-''^ - 2B_ sin^ 0(*) + (3 + cos 26'-''^) {B+ cos2ip + C+ sin 2ip) 
AA cos 6'(*) + S_ (3 + cos 26K*) ) - 2 sin^ 6'(0 {B+ cos 2ip + C+ sin 2(/3) ^ 

4 cos 6l(') (C_ + B+ sin 2(p - C+ cos 2^) 
4A cos 6'(*) + B- (3 + cos 26'(*) ) - 2 sin^ 6'(^) (B+ cos 2^3 + C+ sin 2(p) ^ 

4 cos e'^^ (C+ cos 2^ - B+ sin 2(^ + C_ ) 
4A cos 6'(*) + S_ (3 + cos 26K*) ) - 2 sin^ 6K») {B+ cos 2^9 + C+ sin 2(/3) '' 
(3 + cos 26l(*) ) {B+ cos 2^ + C+ sin 2(^3) - 4A cos 6l(*) - 2B_ sin^ 0^"^ 
4:A cos 6'(*) + S_ (3 + cos 26K») ) - 2 sin^ 6'(') (B+ cos 2(/? + C+ sin 2(^) ^ 



using the definitions: 



^1,2 = 
C'5,6 =(^1,2 
^"7,8 =C'l,2 



I + I/1.2I 



m,2\ 



-1/2 



ni,2 smt^i_2 cos (/? .91,2 - ^1,2 cos 6^^2/1, 2 



Til, 2 smfy^ 2 sm(y9(7i.2 — ^1,2 cos al 2 



(A7a) 
(A7b) 
(A7c) 
(A7d) 

(A8a) 

(A8b) 

(A8c) 

(ASd) 
(A8e) 
(A8f) 



A={CeCj-C5Cs), 
B± = (C2C5 + C2Ce) ± {C2CJ - CiCs), 
C± = {C\Ce + C2C5) ± (C2C7 - CiCs), 

The reflectivity and the emission spectrum and polarization are then determined by eqs. Q, ||SJ), and Q. 

B. COMPLEX ANGLE OF PROPAGATION 

In this Appendix we outline some of the physical properties of the wave propagating in the plasma with complex 
index of refraction (cf. Born & Wolf 1970, §13.2). 
For a medium with complex index of refraction n = riji + inj (where n/j and nj are real), the formal refraction angle 

0*^*-' , as determined by Snell's law, is complex. Let cos f?^*) = (1 — sin^ 5)(*) )i/2 — qqs 9]^ + i cos 9j . Defining the vector 
parallel to the plane of incidence s = (— cos ip, — simp, 0), the wavevector for the transmitted waves can be written: 



kW = !^fsin0«S-cos0Wi 
c V 



c 

UJ 



a(t) 



(*)> 



g(t) 



^(*)^s 



sin 6''^*''s + (ni cos9j — niiCos9ff )z — i{nj cos 9)^ + niiCos9j )z 



(Bl) 



The transmitted electric field has the form E^*^^ oc e^^ ■"" ''^*. Substituting cqn. IBip into this expression, the field 
takes the form: 



E^^^ ocexp 
exp 



g(t) 



(nji cos 9j + nj cos 9j^ )z 



.n(^h 



— ( sin 6''^*'s — {n^cos9)^' — nj cos 9} ')z] — iujt 



(B2) 



Thus, the transmitted wave has a propagating component multiplied by a damping factor. Since the amplitude of the 
wave must decrease as it travels through the medium, eqn. IJB2|) gives the following condition on the index of refraction 
(recall that in the geometry of H3.2I z < 0): 



3(t) 



(«) 



ur cos 9}' + nj cos 9]^' > 0. 



(B3) 
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The traveling component can be used to define a new wavevector k' — sind'-^^s — {tir cos 9 j^ — nj cos 9 j)z. The 

real angle of propagation is then given by cos^*-*' = k' • k'/|k'|. By assumption, the angle of propagation for the 
refracted wave measured with respect to the z axis must be greater than 7r/2, yielding a second condition on the index 
of refraction: 

nW nit) 

ni cos 6^} — 71 _R cos 9}^' 



-KcosO^^y = 



<0. 



sin^ 6l(») + (n/ cos 9y> ~ n^ cos 6*^^ ) 



3(*)^2 



(B4) 



The real and imaginary parts of the indices of refraction for the birefringent transmitted waves must satisfy eqs. (jB3p 
and (El. 
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